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ABSTRACT 


A two-dimensional planer labyrinth seal test rig was designed to operate with air 
supplied at 45 psig and temperatures up to 150 °F. The rig (grates with a manually 
specified test section pressure up to 30 psig yielding Mach numbers to 0.9 and gap 
Reynolds numbers to 100,000. The air flow rate through the seal vdll be controlled by 
setting inlet pressure and adjusting an oudet control valve. The test section measurements 
are 18 inches wide by 1.5 inches depth by 6 inches in length and provides for 10:1 large 
seal" geometry seals to be used to facilitate measurements. Design maximum seal gq> size 
is 0.15 inches. The test section has a glass viewing port to allow flow field measurement 
by non-intrusive means such as Laser Dopplo- Velocimeter (LDV) with seals containing up 
to 5 sealing knives. Measurements of pressure, temperature and flow fields can also be 
simultaneously measured by probes inserted in the seal itself, or mounted on the 
removable/replaceable top plate. Inlet flow is conditioned through the use of a dump 
diffuser incorporating screens, honeycombs, expansion and contraction portions. The inlet 
flow to the test section can be modified from uniform to various non-unifearm conditions by 
employing proflle generators such as screens and winglets. A detailed mechanical design 
has been conducted including stress analysis and seal flow rate predictions. 
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I. INTRODUCTION 


There has been rapid development of increasingly more efficient and power dense 
turbomachinery, particularly in the area of gas turbines. The U.S. Navy is particularly 
intei'ested in improving turbomachinery since most of its active surface combatants use gas 
turbines for main propulsion and electrical power generation. In addition, every aircraft in 
the Navy inventory is driven by some form of gas turbine. Higher efficiency power plants 
directly translate to better range and longer on-station time without refueling, and lower life 
cycle operating costs. Higher power density generally results in smaller machinery and 
increased performance (speed and acceleration). 

Most of the recent advancements in gas turbine efficiency have been realized 
through increased thermodynamic efficiencies of the components, higher cycle pressure 
ratios, and higher turbine inlet temperatures. Signif cant resources are continuing to be 
directed toward improved turbomachinery. For example, a jet engine manufacturer may 
spend tens of millions of dollars and years of research in an attempt to improve the 
efficiency of a gas turbine compressor by only a fraction of a percent. In most cases the 
risk level is very high since there is no guarantee that further development will improve the 
efficiency. While most of the losses in efficiency are due to the losses generated in the 
primary air flow path, some interstage losses occur as a result of ineffective sealing 
between higher and lower pressure sides across each stage. Non-contact labyrinth seals are 
typically employed, and are one method to minimize the interstage losses (and improve 
efficiency) within a gas turbine by reducing this high to low pressure leakage. Figure 1.1 is 
a two dimensional representation of a 3 dimensional annular labyrinth seal. The object of 
this seal is to minimize flow rate by successively throttling the flow in one or more 
chambers to destroy kinetic energy. Most approaches to increasing seal effectiveness 
pertain to closing the sealing gap, adding more "torturous" passages, and increasing the 
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si/e of the chambers. Reduction of interstage losses lead directly to improved compressor 
and turbine overall efficiencies, especially at higher operating pressures. 

Sealing technology is an area now experiencing increased importance. As the 
operating pressures have increased in turbomachinery, their tendency to leak has increased 
since for the same clearance level, a higher pressure ratio leaks more than a lower pressure 
ratio (until the airflow through the seal becomes choked). Some progress has been made in 
non-contact seal technology, but it has significantly lagged behind other technological 
improvements according Wrigley fl|. In the U.S., NASA Lewis Research Center has also 
pursued experimental evaluations to improve sealing technology, the motive for this 
research has been reducing fuel consumption of gas turbine engines for the U.S. Military 
thus reducing operating costs. According to the NASA Lewis report, contained in the 
AGARD collection of 1978, component eflidencies still have potential for improvement, 
and much of that improvement could be achieved through advancements in sealing 
technology. Advanced seals could easily result in a 1/2% efficiency gain for far less cost 
and risk than a complete redesign of the compressor blades (as described above) 11 ]. 

In order to improve sealing technology, a test bed is required to better understand 
the highly complicated flow through a seal, and to adequately collect data in support of 
Computer Fluid Dynamics (CFD) simulations. To provide a facility in support ol 
improving sealing technology, this two dimensicmal frianer pressurized air labyrinth seal 
test rig has been designed. The major goals for this test rig are: 

1. To coiKluct experiments at Reynolds and Mach numbers that are typical of modem 
gas turtxnes on large scale, geometrically similar 2D seals. 

2. Support accurate bulk flow measurements to verify 2D representaticxis of optimum 
non-contact seal designs through modification d* the number of knives, roughness, 
orientation, spacing and geometry. 







3. Make detailed flow Held measurements of sufficient resolution by LDV, hot wire, 
and pressure transducers to compare with CFD code results. 

This repext is organized into five main chapters and is followed by a section of 
riguics containing the test rig design drawings and flow prediction graphs. Cha|Xer I 
provides an introduction, and is fdlowed by Chapter 11-Background, which provides a 
literature review of labyrinth seal research. Chapter Ill-Design Calculations covoa design 
constraints, desired airflow inlet conditions, test rig preliminary design proposals, and 
design of each test rig component including the piping system. Chapter III also contains 
examine calculations. Chapter IV-Data Acquisition System covers the fxoposed computer 
controlled data collection system and its associated components for which the rig is 
designed to emfdoy. Chapter V-Test Rip Performance discusses the factors pertaining to 
labyrinth seal effectiveness and provides a model formulation from which a computer 
program was written to predict Mach number, Reynolds number, and flow rates for 
various seal configurations. This chapter also discusses the results of various simulations. 
Chapter VI-CcxiclusicHis and Recommendatiems describes the overall results of this report 
and provides recommendations for further areas of study. 
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II. BACKGROUND 


The concept of limiting fluid leakage using a number of restrictors is not a new 
idea. In his review of Labyrinth Seal literature, Sneck [2] credits C.A. Parsons with 
development of the labyrinth seal in concert with Parson's [3] development of the steam 
turbine near the end of the 1800's. As Sneck describes it. Parson's idea was to create a 
torturous flow path between high and low pressure regions by means of a series of non¬ 
contacting restrictors and separating chambers. The main function of the restrictors was to 
dissipate the kinetic energy of the fluid flow in the separating chambers and thus reduce the 
leakage from high to low pressure areas. Parson's original method of multiple sealing 
knives has changed very little except for more modem geometry and orientations. 

The first pioneering theoretical paper on labyrinth seal design was introduced by 
Martin [4] in 1908. Martin considered the labyrinth to be a series of discrete throttling 
processes very similar to flow through orifices. Martin derived a formula for labyrinth seal 
leakage based upon a number of simplifying assumptions namely: flow through the seal 
was isothermal, all knives and cavities were symmetrical, kinetic energy was completely 
destroyed in each cavity, and the airflow was sub critical (not sonic, i.e., flow not choked). 
Most researchers following Martin's published work focused on modifying some or most 
of Martin's original assumptions. Gercke [5] modified Marlin’s equation with the addition 
of a kinetic energy carry-over factor. He assumed the flow through the seal was adiabatic in 
that the flow would return to an isothermal condition via a constant pressure process at 
each throttling. This theory supported the establishment of his kinetic energy carry-over 
factor. Gercke also modified Martin's formula for varying areas between the restrictors. 

The next major modifier of Martin's formula was Egli [6] in 1935. Egli examined 
the use of Marlin's formula for fewer than 4 restrictors in both incompressible and semi- 
compressible flows. One of Egli's main contributions was that he noted the pressure ratio 
across each restrictor increases as the flow moves from the first restrictor to tiic last. 







Therefore the last restrictor is always the first (assuming same gap size through the seal) to 
choke (reach cTitical flow pressure ratio). Egli also experimentally determined that the last 
restrictor would not choke at the same pressure ratio as a sharp edged orifice since its flow 
coefficient increases as the pressure ratio increases. Egli produced test results for seals 
showing the effect of gap size and knife thickness on seal performance. 

Hodkinson f7| in 1939 was the first to approach the kinetic energy carry-over factor 
of Gercke from a fluid dynamic standpoint rather than a thermodynamic standpoint He 
considered the interaction of the flow encountering the next sealing knife as a jet with 
portions being sheared off by the obstruction of the knife. Jerie [8] in 1948 examined the 
effect of the knife thickness to clearance ratio on the seal performance. He determined that 
when this ratio was greater than 2:1 than the knife behaves similariy to a rounded nozzle 
instead of a sharp edged orifice. Jerie also determined that a tooth depth to tooth spacing 
ratio of slightly less than 1 yields an optimum performing seal (minimal leakage). There 
were some other papers between 1948 and 1961 of minor importance from Kearton and 
Keh (9) and Zabriske and Sterlicht [10| concerning flow coefficients and friction factors, 
respectively. 

The next major paper was by Vermes (11 ] in 1961. Vermes crmibines the efforts of 
Martin, Gercke, Egli, and Zabriske into one major work. Vermes uses Martin's formula 
and derives new theoretical and semi-theoretical formulas fcx* computation of the leakage. 
These formulas correlate within 5% to experimental tests for three types of seals (straight, 
stepped and comlxnation). Vermes also calculates off (tesign seal performance of the seals 
from both a theoretical and experimental standpoint In 1977 Stocker (12], working for 
NASA Lewis Research Center, investigated the influence of rotation on seal leakage with 
solid-smooth, abradeable, and honeycomb lands. He also made experimental 
measurements of flow leakage for the 2D and 3D cases. His conclusions were that 
advanced seal designs increase seal effectiveness by over 25% when compared to 
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conventional stepped seals, and rotational effects produced mixed results with regard to 
seal leakage depending upon the materials used in constructing the seal. 

The AGARD conference of 1978 [1] specifically states that meaningful work in the 
area of sealing technology has not kept pace with advances in the major gas turbine 
component assemblies, and that resources need to be dedicated to further woric in this area. 

Chapter V-Test Rig Performance fully develops the equations for determining the 
flow through a seal and describes the factors affecting seal leakage, therefore these 
equations will not be developed here. 
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III. DESIGN CALCULATIONS 


A. CHAPTER OVERVIEW 

The test rig design began with first establishing a set of design goals. These goals 
were used to consider a set of design options and tradeoffs and to guide the establishment 
of a set of preliminary designs. The preliminary designs led to selection of one "best 
candidate" design for which a more in-depth stress analysis and cost analysis were 
performed. After ensuring the design had some promise, a complete stress analysis was 
conducted and detailed mechanical design produced. Meanwhile, a set of performance 
predictions was generated to verify the design would achieve the original design set points 
with a low level of risk (both technical risk and safety). 

The sections included in this chapter am Design Goals, Desired Inlet Conditions, 
Design Methodology, Possible Design Approaches, Air supply. Piping and Valves, and 
design of the Diffuser, C^lontraction, Test and Exit sections. The chapter is completed by a 
discussion of the labyrinth seal, and some design configurations used to modify seal 
perfcamance. 

B. DESIGN GOALS 

The first choice to be made was the speciHc design objectives, i.e., what type of 
seals did we desire to test (straight, stepped, planer, annular, etc.) and to what extent did 
we want to test them. In most seals the ratio of gap clearance to rotor diameter is so small 
the 3D annular seal may be considered 2D planer (Figure 1.1), especially if the rotor 
rotation is not particularly high and the fluid friction coefficient is low (as is the case with 
air). With the desire to match gap Reynolds numbers and Mach numbers that occur in real 
gas turbines (with realistic seal geometry and large scale to facilitate instrumentation) it was 
decided that an effective 2D planer seal test rig would suffice. The typical Reynolds number 
range in a gas turbine engine will run to about 10,000 in the compressor and to greater than 
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100,000 in the turbine, with axial Mach numbers less than 0.9 and pressure ratios of 2 to 

3. A 2D annular design (fixed) was considered due to the lack of end wall effects (end wall 
effects are particular to 2D test rigs due to boundary layer growth along the edges), 
however the ability to construct a curved view port was in doubt as was the ease of 
instrumenting the rig (the rig would have to be very large to facilitate instrumentation), this 
rig would be cost |xohibitive. The following majex design goals and are listed below; 

1. Safety was the overriding concern on this project Thus a minimum Factor of 
safety (FOS) of 3 against yield was the design point for each component. 

2. A test rig of the smallest size while preserving good aspect ratio and adequate LDV 
control volume was desired, this resulted in seals scaled upward by about a 10:1 
margin. 

3. Continuous run operation, the airflow supply was limited by available equipment to 
2000 SCFM and 150 psig as shown in Figure 3.1. Blow down operation was 
considered, however this would impact other users of the air supply system by 
depleting the reserve air supply and would require recharge intervals (however 
users of the blow down rig would not significantly affect this rig's operation as it 
would require no recharge time). 

4. Cost and ease of construction were very important, the rig needed to be cost 
effective and "consUaction friendly". 

5. Once constructed, the test rig required maximum expandability, user friendliness, 
and low maintenance. This necessitated the development of a modular design. 
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C. DESIRED INLET CONDITIONS 


Gas turbine engine labyrinth seals, although performing the same sealing function 
in both the compressor and turbine sections, often encounter widely variable (xmditicms of 
airflow conditions based upon design, operating speed, operating temperature and 
pressure. In fact, even seals within the same section such as two sets of comfniessor seals, 
may experience variable inlet air profiles. With this in mind, it was desired for the test rig to 
possess some method able to produce an adjustable air entrance profile at the test secticm 
thus replicating real engine conditions. The efTect of air entrance profiles on seal 
performance is not well understood, thus it is an important aspect to consider in 
constructing labyrinth seals. 

The rig was thus required to produce flow from smooth and uniform prollles to 
turbulent linear and non-uniform profiles, with provisions for introducing shear or vortex 
conditions. A method was also desired to vary pressure ratios from 1.2 to 3.0. 

D. METHODOLOGY 

To provide the range of profiles described in section C. above, a source of airflow 
modification methods was required. Although many books on wind tunnels were 
reviewed, only one particular wind tunnel book by Pope 113] was consulted and its impact 
on the design of the test rig appears in following sections. The wind tunnel information 
was utilized to determine the impact of screens, honeycombs, expansions and contractiems 
on modifying the airflow. Using the information from the wind tunnel literature and 
considering the constraints of part B, 3 designs were considered. These designs were 
evaluated for structural integrity using simple plate equations for stress and strain and 
modified as required. A final design was chosen based upem strength, size, ease of 
construction, and case of modification. Each of the preliminary designs is discussed in 
.section E. 
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E. POSSIBLE APPROACHES 

Proposal #1 consisted of a diffuser designed for maximum pressure recovery and 
included both horizontal and vertical transitions of S-6 degrees with internal splitter vanes. 
The diffuser was designed in 2 sections, a first diffuser stage followed by a fine screen, 
then followed by a second diffuser stage ending with multiple fine screens and a 
honeyoxnb section. The diffuser was expected to require 35 psig in diffuser assemUy and 
employed a 3D large area ratio (6.7) final contraction section would form the air into a very 
smooth profile (less than 2% total turbulence). The problem with this deagn was its 
complexity, projected construction time, size, cost, and possible problems with separated 
or highly ncxi-uniform flow. 

Proposal #2 was similar to #1 with the exception of a modified diffuser assembly 
that would now be one piece with separate 2D horizontal and vertical transitions, an 8 
degree diffusion angle with many s|4itter vanes, and a series of diffuser stage screens still 
ending with 1 honeycomb section. At the diffuser exit would be a smaller contraction ratio 
(3.3 vice 6.7) operate with 35 psig in the diffuser assembly (still designed for maximum 
pressure recovery) and maintain a 3D final contraction section. The major problems with 
this design was complex placement the numerous vanes, still time consuming to 
construct, and most importantly the questionable ability to “shear” the gauss distribution of 
the air inlet into equal parts for effective diffusion. This design would have even more 
chance of being affected by separated ex' highly non-uniform flow than proposal #1. 

Proposal #3 incorporated a diffuser being of a "dump" type (not designed for 
maximum pressure recovery) with a removable top plate for easy repositioning of screeas 
and honeycombs according to guidance received by Dr. Homung [14]. The diffuser would 
contain only vertical transitions with no vanes, but a series of torturous air paths, multiple 
diffuser stage screens, and multifile honeyccxnb sections. Diffuser operating pressure 
would now be 45 psig due to the increased diffuser losses, this was an increase from 35 
psig. The final contraction section changed to only a vertical contraction (Area ratio of 6). 
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The problems with this design are: it is now very heavy, it must be corrosion protected, 
aiKl there is no way to correctly predict the airflow prc^ile and initial perforation plate 
placement. The test section did not chan^ appreciably from the initial deagn, it is 
constructed of steel and 1" glass plate. Sensors are mounted either on seal or cm top plate 
(replacing glass), and are discussed more fully in Chapter IV. 

F. AIR SUPPLY, PIPING AND VALVES 

Air is supplied to the test rig area frcxn a 3 stage 600 Hp Elliot centrifugal air 
compressor producing 2000 SCFM at 150 psig and 180 <>F maximum outlet temperature 
(Figure 3.1). The air passes through an air drier assemUy and to an additional air 
compressor capaUe of boosting the exit pressure to 300 psig. In most circumstances the 
boost pump is bypa.ssed and the air is fed directly into an 8000 ft^ storage tank group that 
supplies the test rig area. The air is piped into the test rig area through 4 ittch diameter 
stainless steel piping, the cormection designed to be utilized by this test rig terminates with 
a steel blank flange with a standard 6 bolt pattern of 1/2" bolts. Since the supply manifold 
runs the length of the building and supplies other test tigs, the upstream building shutoff 
valve can not be used as the sde on/off valve for this test rig since it impacts all air users in 
the building. Therefore, a dedicated 4" full shut off valve rated to 300 psig is required. The 
shut off valve can also be used as a throttling valve when the manifold is pressurized to 300 
psig (when the 2X boost pump is in operation) to prevent pressures greater than 150 psig 
in the test rig piping. 

A drawing of the pipe connection flange, required to replace the manifold blank 
flange, is included as Figure 3.3. This flange is designed for minimal flow loss and 
reduces the 4" manifold diameter to a 2" diameter. To determine the flow loss coefficient 
for this flange design, tables from While f 15] were consulted. Thus for a loss coefficient of 
0.005, a required radius/overall diameter =0.20. Knowing the inlet diameter of the pipe is 
1.939", a value of r is determined to be 0.388". The next larger size (using 1/16" 
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increments) is 7/16” that equals 0.437S*. Utilizing the supplied airflow as a guide (2000 
SCFM), air mass flow was determined from the following equation: 


Ptui^'wr (3-1) 

ft 

{.Cn6^){2iXO^) - 2.55— 
fr min sec 

A pipe diameter of 2" was chosen to reduce the cost of die compcments, and 
facilitate structural design (smaller structures are less massive for equivalent strength). In 
the interest of having a smooth pipe wall with minimal upkeep, and simpliflid pipe 
connections, the pipe chosen for this rig was 2" schedule 80 PVC rated to 600 psig. The 
inside diameter of 2" schedule 80 pipe is nominally 1.939”. PVC jnpe is very temperature 
sensitive, and loses strength with temperatures mudi higher than standard conditions. 
Although the com{xessor oudet temperature can approedi 180 the massive size of the 
storage tanks (8000 ft^ made of 1.5” steel {rfate) and long run of piping exterior to the 
building (apfH^oximately 200 feet) allows the air supply manifdd temperature to remain 
below 110 op, especially with the normal Ntonterey temperatures around 65 op. Schedule 
80 PVC pipe must be derated by 33% strength at 110 op (to 400 psig), and 50% (to 300 
psig) at 140 op according to ASTM D-1784. 

The valves required for this test rig all must be of the full port design, i.e., the 
inside air passage diameter must not be less than the pipe internal diameter otherwise the 
test rig air flow will be reduced and the test rig will fail to achieve designed operating 
points. 
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The following valves (or equivalent) should be procured {nrior to initial construction 
in order to account for any changes required in connections prior to construction, the valves 
should be of sturdy construction and rated to 300 psig (shut off and regulator). Fisher 
series 667-ET stainless steel valves would be a good, but expensive choice: 


Name 

Size 

Type 

Notes 

Manual-Shut-off 

valve 

2" pipe thread to 2" 

pipe thread 

Globe valve 

fuUptnt 

ability to reduce inlet 

pressure 300 to 150 

Adjustable-Pressure 

RegulatcH-valve 

2" pipe thread to 2" 

pipe thread 

Regulator 

full port 

regulate 0 to 45 psig, 

tight shut off. 

Pressure-Relief 

valve 

r'pipe thread 

Relief, spring 

adjustable 

ability to set relief 

range 35-50 psig 

Outlet-Ctontrol 

valve 

4" flange or 4" pipe 

Throttling Butterfly 

or Globe, full port 

max pressure 30 

psig 


Table 3.1 


A filter is included in Figure 3.2, this filter is required since no upstream filtration 
capability exists within the air supply system. The use of the filter will be required to 
ensure long life of all fixed sensors, particularly hot wire probes. The filter chosen is a 
Dollinger coalescing model GP-198-0030-020 rated for 180 psig (maximum) and 840- 
1008 SCFM at 150 psig with a pressure drop of 2 psig. This filter is advertised to have a 
99.9% efficiency at removing .3 micron and larger contaminants (including water and oils). 
A differential pressure gauge is included with the filter assembly and indicates when the 
filter media requires changeout. Inlet and outlet connections are 2" NPT. 
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G. DIFFUSER ASSEMBLY 


To facilitate ease of construction, potential modification, and robust design, a 
rectangular dun^ diffuser is onployed. In cmler to adjust the air profile exiting the diffuser, 
a removable bolted top is included (widi lifting padeyes). The diffuser can be modified with 
a number screens, honeycombs, and perforations if desired. The factOTS govoning the 
application of screens, honeycombs, perforations, and contractimi ratios are discussed in 
section H. Diffuser section drawings are contained in Figures 3.4 to 3.10. 

The Diffuser assembly is comprised of a 1" steel bottom plate with 0.75" welded 
side plates and a 0.75" welded entrance plate. These plates have welded flange connections 
to secure the removable 1" top plate. Internally the diffuser holds welded sections of 0.25" 
steel plate of various lengths to act as spacers for the perforation plates, screens, and 
honeycombs. Allowance is made for two side plates and a bottom plate welded together to 
form each spacer, the exit section flange is stepped down to secure the last set of spacer 
plates. The exit of the diffuser section has a bolted connection flange to the contraction 
section. All bolts used in the diffuser and ccmtraction sectitm are 1/2" diameter SAE grade 1 
or above, coarse threaded. The following factOTS of safety (FOS) exist: 

3.3 diffuser front plate (FOS 9.9) with 3X stress concentration 

factor. 

4.42 diffuser top plate 
9.9 diffuser side plates 
5.75 attachment bolts 

In completing stress calculations for the top plate, an expected maximum operating 
pressure in the diffuser was chosen as 45 psig. The diffuser top plate area subject to this 
internal air pressure was an area 18" wide and 33" in length. Using Roark's Handbook for 
Stress and Strain [16] in the case of a rectangular plate with 4 fixed edges (case 8) the 
following formulas apply; 
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R)r maximum stress at the cotter of the long edge: 


0.4974yfc^ 


(3-2) 


For maximum stress at the center of the plate: 

0 . 2472 <?^»* 

-—; 3 — 

For maximum deflection at the center the [date: 

().02»qb* 


Variables: q - 45psig 

E - 30E6psi 
t ■ thickness 


(3-3) 


( 3 ^) 


From these formulas, the maximum stress at 45 psig on the 1 ” top plate is 7246 psi 
at the edge, 3601 psi at the center, with a maximum center deflection of 0.01”. The same 
formulas are utilized for the surrounding [dates to generate stress predictions. To calculate 
b(dt strength, Shigley's Mechanical Engineering Design book (17] was utilized for bolt 
speciFications. Shigley's book contains specific formulas to optimize the use of bolts and 
fasteners based upcm their stiffness, however his method was very time consuming and 
cumbersome. Through previous experience with these formulas it was determined that by 
choosing a significant factor of safety (generally 4 or more) against the proof yield strength 
of the bolts, a very conservative and yet rapid determination cd* the number of fasteners 







could be made. With this in mind, an initial FOS of 5 was chosen, and based upon 
geometric considerations the resulting number of 1/2" bolts (using the lowest strengdi S AE 
grade 1) was determined to be 34, this in turn resulted in the 5.75 FOS. 

Force on top plate=(45 psi)(18")(33")=26,730 Ibf 
1/2" bolt stress areasO.142 in^ 

Grade 1 yield strengdi=32 kpsi 
(0.142 in2)(32kpsi)=4544 Ib^lt 
so for FOS=5, each bolt can only hold 909 Ibf 
therefore # bolts =26730/909=29.4 bolts 

Screens and honeycombs are the primary flow conditioners used to noodify airflow 
in the diffuser, and they may be arranged at various positions within the diffuser assembly 
utilizing 0.25" spacer plates as discussed previously. Although the perforated sections will 
aid in breaking up the Gaussian air flow {uofile exiting the inlet pipe to the diffuser, no 
reliable data exists concerning the best configuration fOT this task. The preferred location of 
the 1st perforation plate is even with the inlet tube cap at approximately 5.0" from datum. 
Datum is the zero point defined as the inlet wall of the diffuser. The perforated inlet pipe 
attempts to double the area of the mean inlet flow and strain in to the sides of the diffuser. 
The 1st perfnation plate has twice the airflow area of the proceeding perforated section and 
tries to maintain the flow towards the edges of the diffuser. It is located at the 9.25" datum 
location. The second perforation plate again doubles the effective flow area and tries to 
form a somewhat uniform profile flow into the honeycomb section. The second perforated 
plate is located at 9.5" from datum. 
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Screens have seen very widespread use both as turbulence reducers (fine wire and 
spacing far upstream of the test section) and as turbulence enhancers (coarse wire cloth just 
preceding the test secticm). Screens typically behave in the relation as set forth by Dryden in 
1929 [18] ; 


U'^ _ 1 

U, (l + it)’^ 


(3-5) 


where (/],and U\^ are the root mean square of the velocity fluctuations downstream and 
upstream of the screen, respectively in x, y, and z directions. The pressure drop factor k, 
and the number of screens n are also variables. As an example problem, if 2 screens are 
used and each screen is comprised of 18 mesh wire of 0.011" diameter wire with a 
k=0.842 at a velocity of 20 feel/second, the turbulence will be reduced by 46.7%. 

A concern when utilizing screens is the drag on the mesh due to the airflow. In the 
above case the drag/ft^ =kq, where q is the mass flow/ft2, the drag is 0.41b/ft2 and for a 
screen of 9"X18" the total drag would be only 1.13 Ibf on the screen. The recommended 
screens for the test rig should be constructed of 20 mesh 0.017" wire diameter (with 
K=1.8) screen located at the 20", 24", and 28" points from diffuser datum. This spacing 
will allow adequate mixing and evening out of the flow between the screens. These screens 
should yield a tiffbulence reduction of 78.7% using Equation 3-5. 

Honeycombs are primarily used to strain the flow and make it close to ID as 
possible by breaking up large scale structures. TTie honeycomb prevents circulation and 
lateral velocity variations from propagating as the flow progresses. The only real rules of 
thumb for the honeycomb, is that the length of the honeycomb should be 5 to 10 times its 
cell diameter, and the octagon configuration typically yields the lowest pressure drops 
compared to circular or square patterns. The recommended honeycomb is 0.25" diameter 
hexagons of aluminum construction of 2.5" in length. The honeycomb leading edge should 
be located at the 13.5" datum location and secured by press fitting into place and placing 
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0.25” spacer plates upstream and downstream of this location. Wall interfaces should be 
sealed by the use of a flexible RTV or Permatex to prevent errant flow around the 
honeycomb. 


H. CONTRACTION SECTION 

The contraction ratio of any test rig, or wind tunnel, is an important item to consider 
for two reasons: 1) the desired contraction ratio will determine the size of the diffuser and 
test sections, and 2) the contraction ratio has a major effect on the longitudinal velocity 
variations within the flow path. As discussed by Pope [13], a properly chosen contraction 
ratio has a major contribution to improving the performance of a wind tunneiytest rig. To 
determine the effect of contracdon ratio some background is required in accordance with 
Pope's discussion. 

Vu=upstream velocity variation 
vd=velocity variation after the contraction 
Vu=upstream mean velocity 
Vd=mcan velocity after the contraction 
n=contraction ratio 

Writing Bernoulli's equation assuming pressure and gravity loss is negligible 

iP(V. + v.)"=|p(V', + v,)^ (3-6) 


v>2v.V.=vj + 2v,V, 


(3-7) 


Dividing by Vl / Vl and neglecting the (v/V)2 term as being very small 


2k. 



(3-8) 
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(3-9) 


V*/V* -1//I* 



Therefore, the variation ratio of upstream velocities are variable through the inverse 
of the square of the contraction ratio. For a contraction ratio of 6, the mean velocity 
variaticMis would be reduced to 1/36 o( their initial value, or yield a 97.2% reducticxi. The 
contraction section profile was chosen from Pope's book on Wind Tunnel Oesign fl3]. 
Pope also recommends the axitraction cone be led by a settling chamber equal to 0.5 of the 
equivalent air chamber diameter, and followed by settling chamber of the same size just 
pritx to the test section. A high quality inside oxitiaction profile graph is included with 
cotxdinates as Figure 3.12, this figure should be used to shape the ccmtiaction curvature 
with a 9.0” distance between top and bottom plates at the inlet and a distance d* l.S” 
between the same plates at the outlet 

I. SETTLING/MODIFICATION SECTION 

This portion d* the test rig allows the air to "rest” and restabilize after contracting in 
accordance with Pope's discussion in [Bf For this aze of airflow area (1.5”X18”) an 
effective length of 2.94” is recommended. The actual settling camber available length is 5” 
with an additional 1” provided prior to the test section view area. This 3” of available space 
can be used to generate specific profiles as desired by the user. Some specific 
recommended flow patterns zuid how they can be generated are: 

1. Shear [xotile - achieved by placing a fine screen covering only a portion of the 

airflow area and leaving the remaining area free from obstructions. 
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2. Vortex pfoTile - using a series of partial qjan win^ets to generate any number of 
wingtip vortices (generating axial vortex cores). 

3. Turbulent profile - achieved through the insertion d* a series of trip wires to create 
vortex shedding and turbulence, or turbulence generated by the injection of air 
normal to the flow path. 

Example drawings of each profile is contained in Appendix C: Row Modifications. 
The above recommendations are only a few of the many different profiles which could be 
achieved not only through modification of flow in the settling section, but by modification 
of the flow pattern upstream in the diffuser and contraction sections. 

J. TEST SECTION 

The test section is comprised of four major compCMients, two symmetric side plates, 
a bottom {^ate, and a top plate each bolted to one another and belted to the settling and exit 
sections. The side plates can be either machined or welded out of O.S* steel plate, along 
with the top and bottom plates. The recommended construction technique for the side plates 
however, is machining out of 1.5” mild steel, with the botttxn and top i^ates being of 
welded construction. The test secticxi allows for designed gaps of 0.15" with seals 6" long 
and 18" wide. The test section is designed to be positive locking with no adjustment of the 
seal once installed, the seal will be bdted into place from underneath the bottom section. 
Three glass top plates allows the use of a traversing LDV system to measure flow fields 
from full front to back of the seal an virtually all area covering the width of the seal. Each 
glass plate is 1 " thick with top dimensiems of 6”X4”, and bottom dimensions erf* 5.5"X4", 
the fremt and rear surfaces of the plate are beveled to support placement of the {rfates in the 
test section top plate. Each glass |:^ate should be {x-otected frexn direct glass to metal contact 
by the insertion of thin sealing gaskets or putty. Removal of the top plate assembly fa- 
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cleaning or access to the interior of the test section will be accomplished by removing the 
twenty-two 0.375” bolts. While the top plate is removed the glass plates should be adjusted 
to give a smooth surface from the front to back of the top plate. Some f(xm of filler (RTV) 
should be used to fill the small gaps between the glass and steel plates. The experimental 
seal is axistructed of steel to ensure robustness of the edges and ccxners. It can be made in 
1 piece, or multiple pieces depending on the desires of the user. The seal is secured by 
drilling six 0.025” holes in the bottom plate and passing bdts through the bottom plate into 
the seal that itself has been drilled and tapped. Additiorud hedes for the passage of 
instrumentation should be avoided without additional stress factor calculations. 
Instrumentation cables should be passed through to the exit section where they should exit 
the rig at one of the 1 ” back pressure holes and be sealed. Details of the test sectiem design 
are contained in Figures 3.14 to 3.16. The fcdlowing factors of safety (FOS) exist: 


9.9 Flanges inlet and exit 

3.8 Bottom pltde 

6.2 Each plate glass 

3.5 Top plate steel 

12.5 Side plates 

6.75 3/8" attachment bolts 


K. EXIT SECTION 

The exit section is designed to minimize back pressure when exit to atmospheric 
conditions is desired, yet allow back pressure to 30 psig as requited for experimentation. 
Details of this design ate included as Figures 3.17 and 3.18. For operation with minimal 
back pressure, 4 additional 1" diameter tapped holes are provided to vent internal pressure 
if required and to allow passage of the seal instrumentation cables. These holes will be 
plugged for normal operations. 
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IV. DATA ACQUISITION SYSTEM 


A. ACQUISITION SYSTEM OVERVIEW 

To effectively collect, process, store and display large amounts of data from 
multiple sensors an automated data acquisition system was required. NFS has recently 
obtained a number of IBM-PC compatible DATASTOR 486-DX computers with 14" 
color monitors containing 16 megabytes of random access memory and 535 megabyte 
hard drives. These computers have full tower cases and six free expansion slots. The 
school also has numerous National Instruments AT MIO-16F-5 Multi-function I/O 
boards, IEEE 488 General Purpose Interface Bus (GPIB) boards, and SC2070 GPIB 
termination breadboards. The computers are configured with Microsoft Windows V3.1 
and run Lab view for Windows software. Labview is a full featured graphical programming 
system for data acquisition and control, data analysis, and data presentation. Labview is 
easily programmed to accept inputs from plug-in boards and supports numerous data 
collection options beyond the scope of this discussion. This basic system will be interfaced 
with Scanivalve Corporation self-calibrating fixed temperature and pressure transducers to 
provided basic data collection as described in the following sections. This will provide an 
excellent data collection/analysis system. A schematic diagram of the system is included as 
Figure 4.1. 

B. AIRFLOW MEASUREMENT 

To accurately measure the flow through the test rig, a provision for a replaceable 
orifice plate has been made. Orifice units typically contain both upstream and downstream 
pressure taps and utilize a calibrated sonic orifice plates. The plates are replaceable to 
maintain accurate flow measurement through the pressure range of the test rig. Orifice 
plates should generally be chosen to provide less than 1% deviation in the intended realm 
of rig operation pressures. The upstream and downstream pipe connections for the orifice 
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assembly are 2" NPT of 1.939" inner diameter on this test rig. The orifice assembly is 
designed to be mounted just downstream of the pressure regulator valve, and just prior to 
the entrance of the diffuser assembly. In this position the orifice assembly should 
experience pressure less than 50 psig; however, a potential drawback to this location is the 
high pressure drop induced in the orifice which may prevent the test section from achieving 
the desired pressure and flow rate. An alternative position is to mount the orifice assembly 
upstream of the pressure regulator (150 psig side), in this position the orifice may be 
slighdy less accurate due to the larger operating pressures, but will guarantee sufficient air 
inlet to the pressure regulator and test section. The diffuser assembly and pressure regulator 
threaded connections can therefore be easily switched to allow placement of the orifice 
assembly at the location desired while preserving sufficient upstream and downstream pipe 
lengths to yield good orifice results. A straight pipe length of at least 39" should precede 
the orifice assembly for a fully developed turbulent air profile to occur, and a trailing 
straight pipe length of at least 12" downstream of the orifice. Channels lA&B through 
4A&B for both temperature and pressure are dedicated to collecting this data. 

C. FIXED SENSORS 

To determine the condition of the air at any point in the test rig, a series of pressure 
taps and thermocouples and a traversing hot wire are employed. Two computer expansion 
slots are filled with SDIU MK5 digital cards, with each card connected to its associated 
scanivalve. Although the scanivalx ss are driven by the same drive system, one 48 channel 
valve set-up is configured for pressure measurements (channel A), while the other is 
configured for temperature measurements (channel B). Channels 5A through 8B are 
dedicated to the diffuser section, channels 47A&B and 48A&B to the exit section, and the 
remaining 38A&B channels to the test section. A major design deficiency of the test rig is 
its current inability to support rapid repositioning of the sensors while the unit is in 
operation (however, most test rigs do not have this ability). In fact, no repositioning is 
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currently possible without partial disassembly of the test rig. Sensors can be mounted at the 
discretion of the user depending upon the goal of the experiment, but a recommended 
layout is shown in Figure 4.1. In the recommended layout, the left side of the test rig 
contains more in-line sensors than the right side, this will allow the user to determine if any 
ill effects occur due to this placement. Every data location is also supported by at least one 
and sometimes two redundant sensors. These redundant sensors should also ensure 
adequate safety monitoring capability. The scanivalve system is self calibrating and has an 
advertised accuracy of less than 0.1% of full scale. 

A DANTEC traversing hot wire system is also recommended to aid in determining 
flow fields, and as a supplemental mass flow indicator. This system requires two computer 
slots, one for a GPIB timing board, the other for an AT MIO-16F Analog to Digital board. 
A calibration unit is provided to facilitate easy verification of data. This system is software 
controlled and is fully position able through remote means in three dimensions. The results 
of using the hot wire should be compared with the LDV system described in part C. 
Accuracy is advertised as less than 0.1% of full scale. 

D. NON-INTRUSIVE LDV MEASUREMENT 

Many methods exist for determining the state of a flow Held, three of them are hot 
wires, pitot probes/pressure taps, and Laser Doppler Velocimeters. Of the three only one, a 
traversing LDV, can generate high accuracy data exterior to the test rig. The NFS's recent 
acquisition of a dedicated LDV system complete with seeding atomizer promises to yield 
signiHcant results when used in concert with the NFS high speed and low speed test rigs. 
This LDV system is a DANTEC system with an Enhanced Burst Spectral Analyzer, 3 
watt output wtaer cooled laser, and a stand alone data acquisition system. Software is 
provided by the manufacturer to collect, display, analyze, and modify the LDV data. 

The LDV works on the principle of Doppler shifting to determine the ID velocity 
of the flow field. The laser beam is produced, then separated into multiple colored beams 
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(typically split by a prism with the two strongest beams of blue and green collected at the 
output). Each of these beams are then split again, with 1/2 of the beam frequency shifted 
slightly out of phase. Each beam is then focused into the control volume and scattered by 
the motion of the flow field, sometimes seeding material is required to aid in scattering of 
the beams. The beams are then collected in either a forward scatter or back scatter mode 
and analyzed to determine the resulting frequency shifts. Forward scatter is typically the 
better mode when possible, but the test rig will be limited to back scatter operation only. 
These frequency shifts are then converted into velocities and collected to obtain a velocity 
flow field. 
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V. PREDICTED TEST RIG PERFORMANCE 


A. EQUATION DEVELOPMENT 

To accurately predict the leakage rates for a given labyrinth seal geometry, and to 
predict Mach numbers and Reyncdds numbers, a method was required to ^ecdvely predict 
flow rates over large variations in seal gecxn^ry and inlet pressures. In order to accomplish 
this, a program was written in MATLAB (appendix A) to satisfy these requirements. 
Section B discusses reduction of the data into Madi and Reynolds numbers, and the 
speciHc program predictions. In order to develop the flow prediction program certain 
assumptions and simi^ifications were required: 

1. All inlet and outlet conditions of the air are constant, known quantities. 

2. The working fluid (air) is a peifect gas. 

3. Row through the seal is adiabatic. 

4. The air flow behaves in a semi-compressible fashion at speeds below Mach 0.9. 

(i.e. the calculated air density is the average of the inlet and exit air densities across 
the seal) 

5. All jet kinetic energy is dissipated at each intermediate seal cavity. 

6. The kinetic energy carryover factor (P) and vena-contracta modification (Cc) are 
constant for each knife (p=CcP-0.65). 

And although not necessary, the program is written with two fixed parameters: 
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I. The spacing between the sealing knives are constant 


2. The gap clearance of each knife is constant 


The simultaneous flow equations for both ideal incompressible and ideal soni- 
compressible flow will be developed for one sealing knife. Utilizing the mass flow 
equation for ID steady flow: 

m^pAV (5-1) 

and the equation for the density of air 


P- 


P 

RT 


(5-2) 


and the Bernoulli equation neglecting changes in height with work terms set to 7jen\ and 
treating the inlet velocity is zero due to a sufficiently large reservoir (denoting conditions 
upstream of the tirst knife by the subscript ”u” and conditions at the downstream exit by 
"d"): 

(5-3) 


By rewriting terms Equation (5-3) becomes: 


^2(P.-Pj) 


(54) 


Combining equations (54) and (5-2) into Equation (5-1) an equation for the mass flow 
through the seal is obtained: 






To calculate the semi-ixunpressibie How across one sealing knife some 
modifications to Equations (5-1) through (5-5) will result Utilizing the mass flow equation 
(where the represents an average of the upstream and downstream conditions); 


m*^p*AV* 


(5^) 


and the equation for the density of air (assuming the air density is now based on the 
average pressure across the seal): 



(5-7) 


By rewriting terms ot the Bernoulli equation using the values while still denoting 
conditions upstream of the first knife by the subscript "u" and conditions at the downstream 
exit by "d" Equation (5-3) becomes: 



Combining Equations (5-6) and (5-7) into Equation (5-8) an equation for the semi- 
compressible mass flow through the seal is obtained; 


m* - A 


kP^-Pd)(Pu+pj) A 

I RT " 4Wf 


{Pl-Py 


The ideal equation (5-9) must now be modified by a flow coefficient p, which is the 
actual flow through the seal divided by the ideal mass flow. 
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An expresaon for (a is the product the contraction coefficient Cc (which is the 
ratio d* the minimum jet Dow area occurring at the vena-contracta to the minimum area of 
the orifice) and the kinetic energy carry-over factor, p in 117]: 

M^CcP (5-10) 

The theoretical value of Cc is 0.611, but has been experimentally determined to be 0.65 fcH* 
a sharp edged orifice by Lamb ] 19]. An empirical lelatirxi for p based upon seal geometry 
has been determined by Vermes [ 11 ] as: 




1 

(1-a)" 


(5-11) 


8.52 

- 

- + 7.23 
6 


(5-12) 


where I is the distance between seals, and 6 is the gap between the top of the knife and the 
stator with both surfaces static. 

In combining Equation (5-9) with (5-10) a modified flow rate equation can now be 
represented with q as the mass flow per unit width ot the seal. 

« - - py (5-13) 

In reality, the p modification should be 0.65 rex’ the first knife since no kinetic carry-over 
exists, with the product of 0.65 and p fc^ the remaining knives. To simplify equatiem 
development and MATLAB programming, a constant value of 0.65 is used for all values 
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of |i althcxigh the program does calculate and display a "correct* value for p based upon 
Equations (5-11) and (5-12). 

As displayed in Figure 5.1 (a four knife seal). Equation (5-13) can be modified for 
each seal using the upstream and downstream pressures. When setting each q equation 
equal to one-another and canceling the common terms (R,T,p, and & ) the following 
relationship is achieved. 


(Pu2. p,2) = (p,2. P22) = (P22. P32) = (P32. pd2) (5.14) 

Equation (5-14) can now be rewritten into equations for PI, P2, and P3 in terms of only 
Pu and Pd: 



(5-15) 


(5-16) 


(5-17) 


Similar relationships can now be determined for N sealing knives. Since the pressure ratios 
across each knife can now be determined, the use of the compresvsible flow equation for 
isentropic How from White (15] reveals the Mach number across each seal: 





(5-18) 


Using the relationship Ma=V/a, where a=(YRT)-^ (Y=1 -4 for air), the velocity of the air 

How through each knife can be determined. 
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Reynolds numbers can be determined by using the gap 6. the velocity V, and 
kinematic viscosity of air v (with v averaged for the pressure variations across the seal) 
according to the formula 


fi V 


(5-19) 


B. PERFORMANCE PREDICTIONS 

The test rig is designed to operate with inlet pressures as high as 45 psia and outlet 
pressures between 15 psia and 35 psia with seals of 2 to 5 knives and Mach numbers 
below 0.9. The program in Appendix A was validated by several methods for absolute 
pressure ratios to 5 (with atmospheric pressure exit), and for seals with between 2 and 14 
knives. First, the generated mass flow rates at the inlet, outlet and across each knife were 
checked to ensure consistency (they agree to 8 decimal places). Second, intermediate 
pressures were reviewed and Mach numbers checked to ensure increasing Mach numbers 
towards the downstream seals. Third, the program was run and the output compared with 
the curves generated by the well known Martin's seal formula 14] (modified by the flow 
coefficient p-0.65 ), in which the generated curves were virtually exact Fourth, the 
program was run for a g^ clearance of 0.040 inches and 5 knives and the flow rate verses 
pressure ratio [dotted. This graph was compared to a graph previously generated by 
Vermes' (11] in his work comparing a modified version d* Martins formula to 
experimentally derived data Again, the MATLAB program and Martin's fcxmula were 
virtually exact and correlated extremely well to Verme's experimental data (Figure 5.2). 
The MATLAB program was again run for a seal with 14 knives and a pressure ratio cd* 5 to 
determine the intermediate chamber pressures within the seal. These were plotted on the 
graph of pressure ratio verses scaling knife number and compared with a prior graph of 
experimental data (again a product of Vermes' work) fll). The MATLAB program 
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predicted slightly lower chamber pressures than the experimental results and hence a lower 
rate d* seal leakage (which is to be expected due to the assumptions written into the 
program), but the predicted curve differed from the measured curve by only 3% to 5% 
(Figure 5.3). 

Confident that the MATLAB program is reasonably accurate to within 5%, 
Reynolds number and Mach number verses pressure ratio curves were generated for 2, 3, 
4, and 5 knife seals each with gaps of 1/32", 2/32", 3/32", 4/32", 5/32". A series of four 
graphs were produced with exit pressures of 5 psia, 15 psia, 25 psia, and 35 psia (Figures 
5.4, 5.5, 5.6, and 5.7 respectively). Figure 5.8 is actually Figure 5.4 disf^aying only the 
maximum and minimum Reyndds number and Mach number ranges. Figure 5.9 displays 
low speed predicticxis (pressure ratios less than 1.2) for an exit pressure of 15 psia. 

Figure 5.10 shows the effect of varying back pressure on the de-coupling of 
Reyndds numbers and Mach numbers. Thus for a given Mach number, a variation in 
Reyndds number can be obtained by adjustment of the test rig outlet contrd valve to 
produce a back pressure (above atmospheric pressure) on the exit ade of the seal. 
Conversely, if an attachment were connected at the exit to the test section capable of 
drawing a vacuum, the Reyndds number range could be further increased. 

In summary, based on Figure 5.8, the test section should operate with satisfaction 
below Mach 0.9 with Reynolds numbers up to 1(X),000. Figures 5.11 through 5.14 
display the intermediate chamber pressures and seal leakage rates for 23,4, and 5 sealing 
knives as a function d pressure ratio (to 3) and atmospheric exit pressure. Those figures 
graphically show that Mach number is always highest across the last seal since the highest 
pressure ratio exists at that pdnt, independent of the number d* knives utilized. 
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VI. CONCLUSIONS AND RECOMMENDATIONS 


A. CONCLUSIONS 

This theas has presented a design of a test facility to aid in improving sealing 
technology. Based upon the school's available facilities and the level of research desired to 
be ccmducted a list of goals was presented driving the design of the test rig. The 
methodology of design and supporting calculations were presented alcmg with detailed 
mechanical drawings to support construction such a test rig. A proposed data acquiation 
system was presented to allow adequate collection erf* data, and initial operational 
procedures for the rig were included. An overview of the factors impacting labyrinth seal 
performance were presented along with a computer model to predict test rig operating 
parameters. Graphs of the test rig operating realms were produced along with a required 
test rig parts listing. The overall conclusions are: 

1. The 2D test rig design will be cost effective, sim}de to ccmstruct, and provide 
projected performance. 

2. The proposed design met the design goal of operating with a factor safety greater 
than 3.3 at conditions of Regap < 100,000 and Mach < 0.9, to match Reynolds and 
Mach numbers that occur in real gas turbines, and to sup^rt measurements 
significant res(dution for CFD validation. 

4. Modification of the diffuser assembly may be necessary due to the lack of 
information concerning dump-diffusers, and the uncertainty associated with 
perf oration plate, screen, and honeycennb locations. 
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B. RECOMMENDATIONS 

The primary recommendatirxi based upon this thesis is to build a 2D test rig 
according to the provided specifications. Once constructed the test rig should be used to 
validate the MATLAB program flow rate and chamber pressure predictions by experimental 
research, and eventually modify the MATLAB program for variable carry-over factors and 
irregular knife gap and spacing. The rig should then be used to optimize labyrinth seal 
design. After the basic operation of the rig has been substantiated, further modifications of 
the rig should be considered to improve performance as follows; 

1. Consider an eductor system to create a vacuum at the seal exit, this would provide 
additional range for both the available Reynolds number and Mach number without 
raising the internal operating pressure of the system. 

2. Consider perforations and suction on the end walls to minimize boundary layer 
growth and the .subsequent "end wall effect”. 

3. Design and construct a 3D rotational test rig to determine the effects rotation on 
seal performance specifically, this rig could be used for measurements of flow 
coefficients, friction factors, and heat transfer. 

4. Design and construct a modified test section to allow testing of staggered labyrinth 
seals. 
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Figure 3.1 Air Supply Facility Diagram 















3.3 Manifold Connection Flange 
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Figure 3.4 Diffuser Assembly Inlet Plate 
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Figure 3.5 Diffuser Assembly Bottom Plate 
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Figure 3.6 Diffuser Assembly Top Plate 
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Front View , 1" overhang 

33.75" \T 
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Figure 3.8 Diffuser Assembly Right Side Plate 
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Figure 3.9 Diffuser Assembly Perforation Plates 











Side View 

Note: Both plates made of 0.25" 

1st stage perforation plate perforated as drawn. 2nd stage perforation plate 
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Figure 3.10 Diffuser Assembly Perforation Plates 
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Figure 3.11 Contraction Section Inlet and Exit Flange 
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Settling Section Top&Bottom plate Top View 
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Top View 
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Top View 



















































Acquisition System 
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Figure 4.1 Data Acquisition Schematic 
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Figure 5.2 Flow Rate Correlation 












Figure 5.3 Chamber Pressure Correlation 
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Figure 5.4 Test Rig Pcrfonnancc 
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Figure 5.5 Test Rig Performance 
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Figure 5.6 Test Rig Performance 








Test rig performance, Pexlt=35psia 
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Figure 5.7 TCvSt Rig Performance 
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Figure 5.8 Test Rig Performance 
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Figure 5.9 Test Rig Performance 
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Figure 5.10 Test Rig Performance 
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Figure 5.12 Test Rig Performance 
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Figure 5.13 Te.st Rig Performance 
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Figure 5.14 Test Rig Performance 

































































APPENDIX A: SEAL LEAKAGE PREDICTION PROGRAM 

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% 
% This program will determine the 1-0 intermediate pressures and mass flowthrough % 
% a labyrinth seal with 2 knives. The user specifies seal gap, inlet and exit pressures, seal % 
% length and pressure step size. % 

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% 
clear 

% Specify seal gap in inches. 
g=.i5; 

% Specify seal length in inches (Used only for BETA calculation). 

L=2: 

% Specify inlet pressure in psia. 
po=75: 


% Specify outlet pressure in psia. 
pn=i5: 

% Specify pressure step size in psia 
stepsize=2 

N=(po-pn)/stepsize 

% This converts gap and length to meters. 
gm=g*.0254; 

Lm = L*.0254: 

% This converts inlet and outlet pressures to Pascals. 

pom=po’6894.76: 

pnm=pn*6894.76; 

stm = (pom-pnm)/N; 

% This specifies gamma, R and T (in Kelvin) 

R=287; 

T=293: 
mu 1 =.65; 


alpha = 8.52/(7.23 + (Lm/gm)): 
beta=i/((l -alpha)''.5); 
mu=.65*beta; 


fori=1;N, 

pom(i,l) = pnm + (i*stm); 

p1(i,1) = sqrt(.5*pnm''2 +.5*((pom(i,l))''2)); 

q1 (i,1) = ((mul*gm)/sqrt(R*T))*sqrt((pom(i,l))"2-(p1 (i,1))''2); 
q2(i,l) = {(mul*gm)/sqrt(R*T))*sqrt((p1(i,l))''2-pnm''2); 
end 

pom 

pi 

pnm 

qi.q2 
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%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% 
% The following statements modify the previous program for seal knives of 3,4,and 5. % 

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% 

% This is for 3 knives. 

for i=1;N, 

pom{i,1)spnm + (i*stm); 

pi (i,1 )*sqrt((1/3)*pnm''2 + (2/3)*(pom(i, 1 ))*2); 
p2(i,1 )=sqrt({2/3)*pnm''2 + (1/3)*(pom(i,1))^2); 

q1(i.1) = ({mu1*gm)/sqrt(R*T))*sqrt((pom(i,l))>‘2-(pl (i.1))*2): 
q2(i.1) = ((mul*gm)/sqrt(R*T))*sqrt({p1(i.l))^2-{p2(i.1))^2); 
q3(i,1 ) = ((mul *gm)/sqrt(R*T))*sqrt({p2{i,l ))''2-pnm'‘2); 
end 

% This is for 4 knives, 
for i=1:N, 

pom(i,1) = pnm + (i*stm); 

p1(i.1) = sqrt(.25*pnm''2 + .75*(pom(i.1 ))''2): 
p2{i,1 )=sqrt(.5*pnm'‘2 + .5*(pom(i,l ))'‘2); 
p3(i, 1 )»sqrt(.75*pnm'‘2 + .25*(pom{i,1 ))''2): 

q1 (i,1 ) = ((mu1*gm)/sqrt(R*T))*sqrt((pom(i.1))''2-(p1 (i,1 ))''2); 
q2(i.1) = ((mu1*gm)/sqrt(R*T))*sqrt((pl{i,1))^2-{p2(i,1))*2): 
q3(i,1 ) = ((mu1 *gm)/sqrt(R*T))*sqrt((p2(i.1 )K2-(p3(i,1))''2): 
q4(i,1 ) = ((mu1 *gm)/sqrt(R*T))*sqrt({p3(i,1 ))^2-pnm''2); 
end 

% This is for 5 knives, 
for i=1:N, 

pom(i,1) = pnm + (i*stm); 

pi (i.1 )=sqrt(.2*pnm''2 + .8*(pom(i,1 ))''2): 
p2(i,1 )=sqrt(.4*pnm*2 + .6*(pom(i,1 ))''2); 
p3(i,1 )=sqrt(.6*pnm''2 + .4*(pom(i,1 ))''2); 
p4(i,1 )=sqrt(.8*pnm''2 + .2*(pom(i,1 ))''2); 

q1 (i,1 ) = ((mu1 *gm)/sqrt(R*T))*sqrt((pom{i,1 ))''2-(p1 (i,1 ))''2): 
q2(i,1) = ((mu1*gm)/sqrt(R*T))*sqrt((p1 (i,1))''2-(p2(i,1))''2); 
q3(i,1) = ((mu1*gm)/sqrt(R‘T))*sqrt((p2(i,1))'‘2-(p3(i.1))''2): 
q4(i,1 ) = ((mu1 *gm)/sqrt(R*T))*sqrt((p3(i,1))''2-(p4(i,1))''2): 
q5{i.1) = ((mu1*gm)/sqrt(R*T))*sqrt((p4(i.l))''2-pnm''2): 
end 
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Quantity 


20 FEET 


10 FEET 



Quantity 


APPENDIX B: LIST OF MATERIALS 


Airflow piping and valves 


Desciiptiao 


MANUAL SHUT OFF VALVE 2* GLOBE 


FULL PORT 


ADJUSTABLE PRESS REG VALVE 


(M5 PSIG, FULL PORT. TIGHT SHUT OFF 


PRESSURE RELIEF VALVE, ADJUSTABLE 


35-50 PSIG 


EXIT VALVE. 4" BUTTERFLY OR 


GLOBE, FULL PORT 


FILTER ASSEMBLY GP19B0030-020 


Diffuser assembly 


Descripti<» 


1" STEEL PLATE 9.25"X18.5" 


Approximate cost 


(Dcdlais) 


$40 





Approximate cost 


(Dollars) 


$30 



SCREEN/HONEYCOMB ASSEMBLIES 
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Quantity 


Contraction assembly 


DescripdaD 


0.5" STEEL PLATE 9-X17" 


0.5" STEEL PLATE 18-X18" 


0.75" STEEL PLATE 1.5"X22" 


0.75" STEEL PLATE 2.5"X22" 


1" STEEL PLATE 0.75"X 10" 


0.375" X3" BOLTS 


Approximate cost 


(DoUaus) 


$20 


Quantity 


2 


2 


2 


2 


18 


Setting assembly 


Descripdcn 


0.5" STEEL PLATE 2.5"X5" 


0.5" STEEL PLATE 5"X18" 


0.5" STEEL PLATE 5"X19" 


0.5" STEEL PLATE 4.5"X1" 


0375" X3" BOLTS 


Ai^iroxiniate cost 


(Dcdlars) 


$10 


Test Section 


Quantity 


2 


1 


2 I 


_i 

i 
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3 


Description 

Approximate cost 


(Dollars) 

1.5" STEEL PLATE 2.5"X8" 

$30 


0.5" STEEL PLATE 8"X18" 


0.5" STEEL PLATE 7"X 1.5" 


0.5" STEEL PLATE 1"X18" 


0.375" X3" BOLTS 


1"X4"X6" glass plate 


73 

























































74 



























APPENDIX C: FLOW MODIFICATIONS 


SHEAR 
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